GFTeor-UCF 2010-04 

Fermion perturbations in string-theory black holes 

Owen Pavel Fernandez PiedrsQ 
1 Instituto de Fisica, Universidade de Sao Paulo, 
CP 66318, 05315-970, Sao Paulo, Brazil , and 

2 Departamento de Fisica y Quimica, 
Universidad de Cienfuegos, Carretera a Rodas, 
Cuatro Caminos, s/n. Cienfuegos, Cuba 

Jeferson de OliveiraQ 

Instituto de Fisica, Universidade de Sao Paulo, 
CP 66318, 05315-970, Sao Paulo, Brazil 
(Dated: September 13, 2010) 

Abstract 

In this paper we study fermion perturbations in four dimensional black holes of string theory, 
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show that the falloff is identical to those appeared for massless fields in other four dimensional 
black hole spacetimes. 

PACS numbers: 02.30Gp, 03.65ge 



* opavel@ucf.edu. cu 
t jeferson@fma.if.usp.br 



1 



I. INTRODUCTION 



String theory and its generalization in terms of extended objects named M/Dp-branes is 
a promising candidate for a fundamental quantum theory of all interactions |lj , including a 
description of black holes in a quantum gravity framework [2J. 

The scenario of string theory is consistently formulated only in higher dimensions, which 
are usually compactified, so a full description of a black hole whose event horizon is compara- 
ble to the size of extra dimensions has to be in terms of such higher dimensional description. 
So, in order to discuss the production of small black holes in the Large Hadron Collider 
(LHC) and the evaporation of black holes by Hawking effect, we must consider solutions 
coming from a string theory/branes picture. 

One class of solutions which encompass the features of string theory and also can be 
interpreted as black holes are those non-extremal solutions which came from intersections 
of branes in string theory or in M-theory [3 J [4| [5J [6] . The parameters of such black holes can 
be deduced by corresponding compactifications of the higher dimensions. 

In this work we aim at studying the question of stability of a family of four-dimensional 
black holes in string theory. In order to perform this we calculate the quasinormal spectrum 
due to a perturbation of a Dirac field evolving in the background given by the stringy black 
hole. Regarding the stability of extended objects from low energy limit of string theory, 
some work has been done [7], where it was considered a probe scalar field evolving in a 
p-brane geometry and its quasinormal spectrum in order to analyze the stability and the 
relation between the parameters of p-brane and the quasinormal ringing. 

Quasinormal modes (QNM) of black holes has been extensively studied since the pioneer- 
ing work on stability of Schwarzschild singularity performed by Regge and Wheeler [8|. In 
studying quasinormal spectrum, we can gain some valuable information about the param- 
eters which caracterize the solution, because there are definite relations between the QNM 
and the parameters of solution, see [2] and references therein. Such QNM are characterized 
by a well-defined set of complex frequencies and encode the response of black hole to external 
perturbation, so we can study the stability of a black hole against small perturbations due 
to probe fields or the geometry itself trough its quasinormal spectrum. 

Besides, further developments in QNM research in the context of AdS/CFT and holog- 
raphy allow us to calculate the location of the poles of the retarded correlators functions of 
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certain gauge theories [10] and reveals some connection between the dynamics of black hole 
horizons and the hydrodynamics [TT] \1'2\ . 

The paper is organized as follows: Section II gives a review about the casual structure of 
four dimensional black hole obtained from intersections of branes. In Section III we obtain 
the wave equation suitable to analyze the propagation of a Dirac field in the background 
geometry of the four dimensional stringy black hole. Section IV is devoted to numerically 
solve the evolution equation of such fields in the considered background, and in Section V 
we will concerned with the quasinormal stage, and compute the complex quasinormal fre- 
quencies using two methods: Prony fitting of time domain data and the sixth order WKB 
approach. The results of a numerical investigation of the relaxation of Dirac perturbations 
in stringy black holes at very long times are presented in Section VI. In Section VII we 
present an analytical expression for the quasinormal frequencies in the limit of large angu- 
lar momentum, followed by the last section of the paper, which contains some concluding 
remarks of our work. 

II. (3+l)-DIMENSIONAL BLACK HOLE SOLUTION FROM INTERSECTING 
BRANES 

The metric of the non-extremal 4-dimensional black hole obtained from intersections of 
five-branes read as [I] 



This metric represents a family of black hole solutions parametrized by r#, which gives 
the location of the event horizon, and four charges Qi, Q 2 , Q3, Qi given by : 



that are written in terms of five-brane parameters 5i, resulting from the compactification of 
higher dimensions. This solution describes a four dimensional Schwarzschild black hole in 
the case which the charges Qi are all zero. Also the metric is asymptotically flat and, as we 




(1) 



where 





Qi = sinh 2 (5i), i = l,2,3,4 
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shall see below, has a regular event horizon located at r = r# even with all charges different 
from zero. Besides, the redshift goes to infinity in such surface. 

Infinity redshift surfaces can be found, for a given metric, trough the following relation 

90o( X source) /n\ 

V 9oo(x a ) 

which relates the frequency v measured by an observer at rest away from the source, whose 
emission frequency of say, light pulses, is u . 

In order to have a infinity redshift surface, the frequency v has to be zero, it means the 
frequency v was infinitely delayed due to gravitational effects. So, equation^ implies 

#00 source) = 0; ( 4 ) 

given us the location of the infinity redshift surfaces. In the case of the spacetime of four 
dimensional stringy black holes ([T| it is easy to see that r = r# and r = are the surfaces 
where the redshift of in-falling objects goes to infinity. 

As in Schwarzschild solution, such surface r = r# acts as a one-way membrane for physical 
objects whose trajectories lie in or on the forward light cone. To see this, suppose a smooth 
hypersurface S defined by the equation 



u(x a ) = const. 

The vector iV a = d a u is normal to S, so if S is a one-way membrane it has to be a null 
hypersurface, therefore the norm of iV a has to be null as well. So, 

g ab N a N b = 0, 

determine the one-way membranes, which in our case are located at r = th and r = 0. 
However, for the one of the four charges Qi equal to zero, only r = r# is a null and infinity 
redshift hypersurface, r = is a genuine space-time singularity, as we will see below in more 
detail. 

Let us consider the Kretschmann invariant for the spacetime Q, where we take for 
simplicity Qi = q (the result holds as well for four different charges): 

R abcd Ra bcd = t ^i^(r), (5) 

{r H q + r) 8 



where 



P(r) = r A H q 2 {h + Aq + q 2 ) - 2r 3 H qr(3 + 7q + 2q 2 ) + (r H r) 2 (3 + 8q + 8g 2 ) - 2r H r 3 + r 4 . 
As we see from the above expressions, for r — > r# we have 

R abcd R ab cd(r -> r H ) -> const, 

and for r —7- 

# abcd iW^ 0) -> [r|rg 2 (5 + 4g + g 2 )] . (6) 

We see clearly that r = is not a spacetime singularity but, as previously pointed out 
by Horowitz et al [13] it is actually a inner horizon. However, if at least one of the charges 
is zero, the surface r = becomes singular as we see explicitly in the expression (|6]). 

We conclude that for at least one charge Qi zero, we have a four dimensional black hole 
with a spacetime singularity located at r = covered by a event horizon at r = Th which is 
also a infinity redshift surface, as expected for a spherically symmetric space time. On the 
other hand, if all charges are non-zero, we have a regular black hole with a event horizon at 
r = Tn and a inner horizon at r = 0. 

III. FUNDAMENTAL EQUATIONS 

In curved space-time the massless Dirac equation is written as 

= o, (7) 

where y = T^V^ is the Dirac operator that acts on the four-spinor Y ^ are the curved space 
Gamma matrices, and the covariant derivative is defined as = — \u a ^ a lb-, with \i and 
a being tangent and space-time indices respectively, related by the basis of orthonormal one 
forms e a = e° . The associated conection one-forms uj a ^ = u ab obey de a + u a b A e b = 0, and 
the 7° are flat space-time gamma matrices related with curved-space ones by T M = e^ a 7 a . 
They form a Clifford algebra in d dimensions, i e, they satisfy the anti-conmutation relations 
{ 1 a ,l b } = -2r] a \ with r/ 00 = -1. 
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First consider some properties of the Dirac operator that will be used in future calculations 
[TU [T5] . Under a conformal transformation of the metric of the form : 

g^y = fi 2 ^, (8) 

the spinor ip and the Dirac operator transforms as 

^ = fi-§$, (9) 

y^ = fHy^. (io) 

If the line element of a space-time metric takes the form of a sum of independent components 
ds 2 = ds\ + ds 2 ,, where ds\ = g a b(x)dx a dx b and ds\ = g mn (y)dy m dy n then Dirac operator y 
satisfies a direct sum decomposition 

f=f x +fy (11) 

The above mentioned properties of spinors and the Dirac operator in direct sum metrics and 
conformal related ones, allows a simple treatment of the massless Dirac equations in curved 
spacetimes. For two conformally related metrics, the validity of massless Dirac equation in 
one implies the the validity of the same equation in the other. Then, the idea is to solve 
the Dirac equation in the curved space described by an initial metric tensor performing 
successive conformal transformations that isolate the metric components that depend of a 
given variable, and applied successively the direct sum decomposition of Dirac operator, 
until to obtain an equivalent problem in a spacetime of the form M 2 x J^ 2 , where M 2 is 
two-dimensional Minkowski spacetime and J^ 2 is a two-dimensional maximally symmetric 
one, in which the spectrum of massless Dirac operator is known. 

The method work as follows. First, in the following we suppose that our four dimensional 
metric is spherically symmetric, and given by: 

ds 2 = g^dx^dx u = -A(r)dt 2 + B(r)dr 2 + C(r)dfi 2 , (12) 

where dQ 2 denotes the metric for the (2)-sphere S 2 . 

It is easy to show that with the identification (t, r, 8, 0) — > (0, 1, 2, 3) if we choose the basis 
one-forms as e° = e°(t, r) = f^(t,r)dt+ f( r )(t,r)dr, e 1 =e 1 (t,r) = g( t )(t,r)dt+g( r )(t,r)dr 
and e k = e k (8,(p) = h^(8, <p)d8 + h?X{6, 4>)d4>, k = 2, 3 then we have some conection one- 
forms that are equal to zero, and only oj° q , u\, and u/„-; i, j = 2, 3 are different from 
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zero. This fact allows us to write V a = V I(2)j a = 0--3, where 1(2) is the unit matrix in two 



dimensions 

Now, under a conformal reescaling of the form ds 2 = C(r)ds 2 , where: 



ds 2 = 


i 

A 

= ~c 


dt 2 + 


— dr 2 + dtt 2 2 , 


(13) 


we need to solve the equation ^fijj 


= 




~ 3 

= 0, with ^ = C*if). 


In general, the Dirac 


matrices can be choosen as: 












f° = 


" 1 (2) 


® 1(2), 


(14) 




f 1 -- 


" 1 (2) 


® 1(2), 


(15) 




f 2 = 




^ r° 2 ), 


(16) 




f 3 = 


- f 5 
" 1 (2) 


®f(2), 


(17) 



where r? 2 N, rL-, and TjL = — r^TjL are two-dimensional Gamma matrices, and (r|L) 2 



1(2). 



Since the orbit-space part and the angular part of the metric are completely separated, 
one can write the Dirac equation in the form: 



(f ° 2) V + f f 2) V x ) ® I (2) + f f 2) ® (f a V a ) 



s 2 



^ = 0, (18) 



where ^r a V a j denotes de Dirac operator in the 2-sphere, whose orthogonal set of eigen- 



S 2 

spinors |^ are defined by [T7] : 



f a V a ) g ±) = ±i(^ + l)g ±) , (19) 
where / = 0, 1, 2, . . . . Now expanding ip as: 

^ = E(^ +) +^ _) )- (20) 

i 

we can put the Dirac equation in the form: 

{f ? 2) Vo + f J 2) Vi + f%, [±i (J + 1)] } <pf > = 0, (21) 
In the following we shall work with the + sign solution, because the — sign case can be 



,~ 2 



worked in the same form. Denoting as ds\ = — ^dt 2 + ^dr 2 the t — r part of the metric 



G 1 cr 



(13) and doing a conformal rescaling of it in the form: 
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d~ Sl 2 = - [-dt 2 + drl] , 
where dr\ = J 4dr is the tortoise coordinate, we have: 



(22) 



1% + ^d^+t^J^ (£+1) 



= 0, 



(23) 



being 7* and Y* Dirac matrices in the two-dimensional Minkowski space-time with metric 
ds 2 = -dt 2 + drl 

Now choosing the representation of Dirac matrices given by: 



(24) 



where the a 1 are the Pauli matrices: 



a 



1 

1 



-i 

1 



a 



1 
-1 



(25) 



and writting 

Xt(t,r) 

we obtain the following equations for each component of the Dirac spinor iff. 



(26) 



and 



where 



■ d(e ,9X1,. n 
% ~dt + dr~ + ai = ' 



^— - \-A £ Q = 0, 



dt <9r* 



A/(r) 



(£+1). 



The above equations can be separated to obtain: 



d 2 Ce d 2 Q 



and 



where: 



dt 2 dr 2 

d 2 Xi _ d 2 xe 
dt 2 dr 2 



+ V + (r)Ce = 0, 



+ V-(r) X e = 0, 



dA ( 
dr* 



V± = ±^ + A 2 . 



(27) 

(28) 
(29) 

(30) 

(31) 
(32) 
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Figure 1. Potential V(r) for Dirac perturbations with t = (bottom) to I = 5 (top). 

The above equations gives the temporal evolution of Dirac perturbations outside the black 
hole spacetime [18J. As the potentials V + and V_ are supersymmetric to each other in the 
sense considered by Chandrasekhar in [Tj5], then Q(t,r) and Xt{tjf) wu l have similar time 
evolutions and then they will have the same spectra, both for scattering and quasi-normal. 
At this point it should be stressed that for the spinor (p\ , we have these two potentials 



again. In the following we will work with equation (30) and we eliminate the subscript (+) 
for the effective potential, defining V(r) = V + (r) . 

Now for (3+l)-string theory black holes with line element 0, making the identifications 
A(r) = f^ 2 {l- r -f), B(r) = f 1 ' 2 (l - and C(r) = r 2 / 1/2 , we can calculate V{r) 



using (32) with A^ given by: 



k,= f —Jl- — (t+l). (33) 



r 



In Figure ([T| we show the effective potential for massless Dirac perturbations in four di- 
mensional stringy black holes for various multipole numbers I. In this Figure distances are 
measured in units of the black hole horizon radius r#. We see that V(r) has the form of a 
definite positive potential barrier, i.e, it is a well behaved function that goes to zero at spa- 
tial infinity and gets a maximum value at a well defined peak near the event horizon. Then 
we can expect that the stringy black holes are stable under massless Dirac perturbations, a 
fact supported by the numerical results that we will present in the next section. 



IV. TIME EVOLUTION OF DIRAC PERTURBATIONS 



In order to integrate the equation (30) numerically we use the technique developed by 
Gundlach, Price and Pullin [20] ■ The first step is to rewrite the wave-like equation (30) in 



w 

-9- 



N 



C 





-*E 



Figure 2. A cell of the integration grid in the plane (u, v) limited by the points N, S, E, W e C. 
This cell represents an integration step, and the initial data are specified on the left and bottom 
sides of the rhombus. 

terms light-cone coordinates du = dt — dr+ and dv = dt + dr+ 

:— + V(u,v)jUu,v)=0. (34) 

and use the following discretized version of the above equation 

An An 

(e(N) = Q(W) + Q(E) - Q(S) - -^V(S) (Q(W) + Q(E)) + 0{h% (35) 

where the letters S, W, E, N are used to mark the points that limits a particular integration 
cell of the the grid according to: S = (u, v ), W = {u + Au,v), E = {u,v + Av), N = 
(u + Am, v + Av). In figure [2] we show the cell of a given integration step. We see that the 
field value at point iV depends only of the field values at the points S, E and W. Given a 



set of initial conditions at the two null surfaces u = uq and v = vq, we can find, using (35), 



the value of the field Q inside the rhombus which is built on this two null surfaces. Then, 
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by iteration of the integration cell, we find the complete data describing the evolution of the 
fields with time. 

The obtained results from the integration in the case of massless Dirac fields in stringy 
black hole background can be observed as the time-domain profile showed in figures ^ to 
(l7j). In such profiles r = 3r# and the time is measured in units of black hole event horizon. 



0.2 




50 100 150 200 250 20 50 100 200 500 

t/r H t/r H 

Figure 3. Normal (left) and logaritmic (right) plots of the time-domain evolution of £ = massless 
Dirac perturbations. 




' 50 100 150 200 250 10 20 50 100 200 500 

t/r H t/r B 

Figure 4. Normal (left) and logaritmic (right) plots of the time-domain evolution of £ = 1 massless 
Dirac perturbations. 



As we can see, the temporal evolution of Dirac perturbations follows the usual dynamics 
for fields in black holes spacetimes. After a first transient stage strongly dependent on the 
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50 100 150 200 250 20 50 100 200 500 

t/r H t/r H 

Figure 5. Normal (left) and logaritmic (right) plots of the time-domain evolution of £ = 2 massless 
Dirac perturbations. 




Figure 6. Normal (left) and logaritmic (right) plots of the time-domain evolution of £ = 3 massless 
Dirac perturbations. 

initial conditions and the point where the wave profile is observed, we see the characteristic 
exponential damping of the perturbations called quasinormal ringing, followed by a so-called 
power law tails at asymptotically late times. 

It is important to mention that we perform and extensive numerical exploration of the 
perturbative dynamics for different values of angular number £, and we have not found any 
instability for massless Dirac perturbations in four dimensional stringy black holes. It is 
an expected result due to the positive definite character of the effective potential, as we 
mentioned in the previous section. 
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V. QNMS USING SIXTH ORDER WKB METHOD AND PRONY FITTING OF 
CHARACTERISTIC DATA 

In the following we will assume for the function Q(t,r) the time dependence: 

Q(t, r) = Z t (r) exp(-iuet), (36) 

Then, the function Z^{r) satisfy the Schrodinger-type equation: 

d 2 Z 



drl 



+ [u 2 - V(r)] Z e (r) = 0. (37) 



The quasinormal modes are solutions of the wave equation (30) with the specific boundary 
conditions requiring pure out-going waves at spatial infinity and pure in-coming waves on 
the event horizon. Thus no waves come from infinity or the event horizon. 

In order to evaluate the quasinormal modes we used two different methods. The first 



is a semianalytical method to solve equation (37) with the required boundary conditions, 
based in a WKB-type approximation, that can give accurate values of the lowest ( that is 
longer lived ) quasinormal frequencies, and was used in several papers for the determination 
of quasinormal frequencies in a variety of systems [7J [181 I2TH25] . 

The WKB technique was applied up to first order to finding quasinormal modes for the 
first time by Shutz and Will [21 J. Latter this approach was extended to the third order 
beyond the eikonal approximation by Iyer and Will [22] and to the sixth order by Konoplya 
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t 


n 


Sixth order WKB 


Prony 








0.1527- 0.0623i 


0.1525 - 0.0620? 


1 





0,3147-0.0628?' 


0,3147- 0.0626? 


2 





0.4744 - 0.0629i 


0.4744 - 0.0629?' 


2 


1 


0.4663 -0.1899? 


0.4661 - 0.1896? 


3 





0.6337 - 0.0629i 


0.6337 - 0.0629?' 


3 


1 


0.6275 - 0.1894? 


0.6271 - 0.1880? 


3 


2 


0.6156 -0.3180i 





4 





0.7928 - 0.0629i 


0.7928 - 0.0629?' 


4 


1 


0.7878 -0.1892? 


0.7878- 0.1892? 


4 


2 


0.7781 - 0.3169? 


— 


4 


3 


0.7641 - 0.4466i 





5 





0.9518 - 0.0629i 


0.9518- 0.0629? 


5 


1 


0.9476 -0.1891? 


0.9476- 0.1891?' 


5 


2 


0.9394 - 0.3162? 




5 


3 


0.9275 - 0.4448? 




5 


4 


0.9121 - 0.5755?' 





Table I. Dirac quasinormal frequencies loth for t = to i = 5. 



J. We use in our numerical calculation of quasinormal modes this sixth order WKB 
expansion. The sixth order WKB expansion gives a relative error which is about two orders 
less than the third WKB order, and allows us to determine the quasinormal frequencies 
through the formula 

■Sf-r >+ ^ <38) 

where n = 0,1,2,... if Re(cu) > or n = —1,-2,-3,.. if Re(u) < is the overtone 



number. In (38) Vo is the value of the potential at its maximum as a function of the tortoise 
coordinate, and Vq represents the second derivative of the potential with respect to the 
tortoise coordinate at its peak. The correction terms IK depend on the value of the effective 
potential and its derivatives ( up to the 2i-th order) in the maximum, see [27] and references 
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e 


n 


Sixth order WKB 


Prony 


6 





1.1107-0.0629* 


1.1108 - 0.0629* 


6 


1 


1.1071 - 0.1890* 


1.1071 - 0.1890* 


6 


2 


1.100-0.3158i 





6 


3 


1.0897 - 0.4437i 





6 


4 


1.0762 - 0.5732i 


— 


6 


5 


1.0600 - 0.7045i 





7 





1.2696 - 0.0629* 


1.2697- 0.0629* 


7 


1 


1.2665 - 0.1890 


1.2665 - 0.1889* 


7 


2 


1.2603-0.3155* 


— 


7 


3 


1.2511 - 0.4430* 


- 


7 


4 


1.2392 - 0.5716i 




7 


5 


1.2246 - 0.7018* 




7 


6 


1.2077-0.8336* 





Table II. Dirac quasinormal frequencies uirn for £ = 6 and 1 = 7. 



i 


n 


Sixth order WKB 


Prony 


8 





1.4285 - 0.0629* 


1.4287-0.0629* 


8 


1 


1.4257-0.1889* 


1.4258- 0.1888* 


8 


2 


1.4202 - 0.3154*' 




8 


3 


1.4120 - 0.4425* 




8 


4 


1.4012 - 0.5706* 




8 


5 


1.3881 - 0.6999* 




8 


6 


1.3727- 0.8305* 




8 


7 


1.3554 - 0.9628* 





Table III. Dirac quasinormal frequencies urn for i = 8. 



therein. 

The second method that we used to find the quasinormal frequencies was the Prony 



t 


n 


Sixth order WKB 


Prony 


9 





1.5873 - 0.0629* 


1.5876- 0.0629* 


9 


1 


1.5848- 0.1889* 


1.5850- 0.1888* 


9 


2 


1.5798-0.3152i 


— 


9 


3 


1.5724 - 0.4421* 





9 


4 


1.5627-0.5698* 


— 


9 


5 


1.5507- 0.6985* 


— 


9 


6 


1.5367- 0.8283* 




9 


7 


1.5207- 0.9594* 




9 


8 


1.5030- 1.0920* 





Table IV. Dirac quasinormal frequencies corn for I = 9. 



i 


n 


Sixth order WKB 


Prony 


10 





1.7462 - 0.0629* 


1.7466- 0.0629* 


10 


1 


1.7439- 0.1889* 


1.7442 - 0.1887* 


10 


2 


1.7393- 0.3151* 




10 


3 


1.7326 - 0.4419* 




10 


4 


1.7237- 0.5692* 




10 


5 


1.7127- 0.6974* 




10 


6 


1.6998 - 0.8266* 




10 


7 


1.6850 - 0.9569* 




10 


8 


1.6686- 1.0884* 




10 


9 


1.6505 - 1.2212* 





Table V. Dirac quasinormal frequencies urn for £ = 10 . 



method [22, E7] for fitting the time domain profile data by superposition of damping expo- 
nents in the form 

p 

k=l 

Assuming that the quasinormal ringing stage begins at t — and ends at t — Nh, where 



0.0 


1= 1 


2 3 


4 5 


6 7 8 9 10 


12 3 


4 5 


6 


7 


8 9 10 




n = - 1 ■ • 


V A 


□ 


O T A ♦ ■ 


• ■ ♦ A 


T 


u 




A V . n _ 


-0.2 


n = -2 ■ • 


V A 


o a 


OTA 


♦ A 


T 


□ 





a v • n = 1 




« = — 3 ■ • 


V A 


<> □ 


O T 




▼ 


□ 





A V • (J = 2 


-0.4 


n = -4 ■ • 


V A 


□ 


o 




▼ O 






a v • n = 3 


-0.6 


n = — 5 ■ « 


V A 















a v • n = 4 




n = -6 • 


• V A 










□ 


• 


A v • « = 5 


-0.8 


n = -7 ■ 


• V A 














A v • H = 6 


-1.0 


« = -8 ■ 
« = -9 ■ 
















A v % n = 1 
v • n = 8 


-12 


n = -10 ■ 
















. re = 9 























-2-10 1 2 

Re u» 



Figure 8. Massles Dirac quasinormal frequencies for various values oil e n. 



N > 2p — 1, then the expession (39) is satisfied for each value in the time profile data 

p p 
x n = i, (nh) = C k e~ lWknh = C k z n k . (40) 

k=l k=l 

From the above expression, we can determine, as we know h, the quasinormal frequencies 
bJi once we have determined Z{ as functions of x n . The Prony method allows to find the 
as roots of the polinomial function A(z) defined as 

p p 

A(z) = \[(z-z k ) = Y j amZ p - m , « = 1. (41) 

fc=l m=0 

It is possible to show that the unknown coefficients a m of the polinomial function A(z) 
satisfy 

p 



m=l 



Solving the N — p + 1 > p linear equations (42) for a m we can determine numerically the 
roots z a and then the quasinormal frequencies. 

It is important to mention the fact that with the Prony method we can obtain very 
accurate results for the quasinormal frequencies, but the practical application of the method 
is limited because we need to know with precision the duration of the quasinormal ringing 
epoch. As this stage is not a precisely defined time interval, in practice, it is difficult to 
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determine when the quasinormal ringing begins. Therefore, we are able to calculate with 
high accuracy only two or, sometimes three dominant frequencies. 

Tables (I) to (V) shows the results for the quasinormal frequencies measured in units of 
black hole horizon radius r# for some values of the multipole moment i. As it is observed, 
the sixth order WKB approach gives results in good correspondence with that obtained by 
fitting the numerical integration data using Prony technique. 

As usual, the oscillation frequency increase for higher multipole and fixed overtone num- 
bers. The fundamental mode, i.e, that with i = and n = is more long-lived with respect 
to the other modes, but it is interesting that imaginary part of the n = quasinormal 
frequencies reach very quickly a fixed value for higher multipole numbers. 

This situation is different for higher overtone numbers, and the imaginary part for a fixed 
n increases with £. For n = 1 we see that its reach a constant value from the i = 8, within 
numerical accuracy. As we expected for stability, all quasinormal frequencies calculated in 
this work have a well defined negative imaginary part. 



VI. LONG TIME TAILS 

Another important point to study is the relaxation of the perturbing fermion field outside 
the black hole. It is a known result that in Schwarzschild black hole neutral massless fields 
had a late-time behavior for a fixed r dominated by a factor t~( 2 ^+ 3 ) for each multipole 
moment I [281129]. 

To study the late-time behavior, we numerically fit the profile data obtained in the 
appropriate region of the time domain, to extract the power law exponents that describe the 
relaxation. As a test of our numerical fitting scheme, we obtained the power law exponents 
for the massless Dirac field considered in this paper in the space-time corresponding to four 
dimensional Schwarzschild black hole with unit event horizon. As we expected, the results 
obtained are consistent with the power law falloff mentioned in the previous paragraph. 

The case for the late-time fermion decay in (3+l)-stringy black hole is similar, as shows 
the results of our numerical fitting presented in Figures @ to tg. 



For I = modes, the numerical decay factor is proportional to t 3,09 , for I = 1 is 
proportional to t~ 5 ' 03 . From the results displayed in Figures @ to (g, and the other 



obtained fitting the time domain data for all the values of angular numbers considered 
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150 200 300 500 700 

t/l-K 

Figure 9. Tail for £ = 0. The power-law coefficients were estimated from numerical data represented 
in the dotted line. The full red line is the possible analytical result. 
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Figure 11. Tail for 1 = 2. The power-law coefficients were estimated from numerical data repre 
sented in the dotted line. The full red line is the possible analytical result . 
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Figure 12. Tail for 1 = 3. The power-law coefficients were estimated from numerical data repre 
sented in the dotted line. The full red line is the possible analytical result . 
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in our numerical work, we can determine the power law factor that dominates the falloff 
of the Dirac field at very late times, in the space-time of the stringy black hole. If we 
suppose that the decay in the Dirac case is governed by factors of the form t-O^+z 3 ) for 
each multipole moment £, then using our numerical results we obtain a decay factor for 
the Dirac perturbations at late times, that is of the form t~( 2e+3 \ We remark at this point 
that this dependance is only a result consistent with our numerical data for all values of 
the multipole moments studied, and remains an open problem to obtain this directly from 
analytical calculations. Then, we can conclude that, outside four dimensional stringy black 
holes, as well as Schwarzschild black hole, the massless Dirac field shows identical decay at 
late times. 



VII. LARGE ANGULAR MOMENTUM 

Although in general the calculation of quasinormal frequencies is performed though nu- 
merical tools, we can use the WKB formula to get an analytic expression for the frequencies 
in the limit of large angular momentum £ . It can be obtained expanding the effective po- 
tential in terms of small values of \ji and using the lowest order of WKB approximation. 
In the following we suppose for simplicity Qi — Q. So, we get 



lx / „d 2 a(r. 

w' = fcr{r m ) — i{n + 

where 



W 2 = fV(r m ) 1)^-2^=, (43) 



and r m is where the peak of V(r) occurs, and its given by 




3^ 2 



The quasinormal frequencies calculated using the above expression completely agree with 
those ones obtained through the six-order WKB approach in the large angular momentum 
regime. In the table (VI), we can see that the imaginary part of the frequencies for fixed 
overtone number n = goes to a constant value in the limit of large values of £. On the 
other hand the real part increases as £ becomes larger. 
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1 


Dirac quasinormal frequencies 


50 


3.96989 - 0.0314673i 


100 


7.93959 - 0.031468H 


150 


11.9093-0.0314682* 


200 


15.8791 - 0.0314683i 



Table VI. Dirac quasinormal frequencies in the large angular momentum limit with n = 0, ru = 2 
and Q = 1. 

VIII. CONCLUDING REMARKS 

We have studied the evolution of massless Dirac perturbations in the space-time of a 
(3+l)-dimensional black hole solution from string theory. Solving numerically the time 
evolution equation for this perturbations, we find similar time domain profiles as in the case 
of Dirac fields in other four dimensional black hole backgrounds. At intermediary times the 
evolution of fermion perturbations is dominated by quasinormal ringing. We determined 
the quasinormal frequencies by two different approaches, 6th order WKB and time domain 
integration with Pronny fitting of the numerical data, obtaining by both methods very close 
numerical results. 

At very late times, the evolution of fermion perturbations in four dimensional stringy 
black holes shows a power law falloff proportional to t~( 2e+3 \ as functions of the multipole 
number. This behavior is identical to that appeared in the late time evolution of Dirac fields 
in other four dimensional black holes. 

There are extensions of this work that would be considered in the future. In the first place 
it would be interesting to study of Dirac perturbations in higher dimensional stringy black 
holes. Another interesting problem is related with the analytical investigation of the late- 
time behavior of massless Dirac perturbations outside stringy black holes, to obtain a result 
in correspondence with our numerical estimate for the decay factor. We also find interesting 
the consideration of the influence of vacuum polarization effects in the quasinormal spectrum 
for semiclassical stringy black holes. The solutions to some of the above problems will be 
presented in future reports. 
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